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Abstract 

We study non-Gaussianity induced by a pseudo Nambu-Goldstone boson with a 
cosine-type scalar potential. We focus on how the non-Gaussianity is affected when 
the pseudo Nambu-Goldstone boson rolls down from near the top of the scalar 
potential where the deviation from a quadratic potential is large. We find that the 
resultant non-Gaussianity is similar to that obtained in the quadratic potential, if 
the pseudo Nambu-Goldstone boson accounts for the curvature perturbation; the 
non-Gaussianity is enhanced, otherwise. 



1 Introduction 



Currently, non-Gaussianity of the cosmological perturbation is a hot topic, since Yadav 
and Wandelt reported a detection of large non-Gaussianity in the WMAP Syr data on the 
cosmic microwave background (CMB) temperature anisotropy [Ij. The latest WMAP Syr 
data is consistent with vanishing non-Gaussianity at 95% C.L. pj, although the data still 
seems to favor similar deviation from the pure Gaussian statistics. It is not yet settled 
whether there are non-Gaussianities, in spite of many recent activities searching for their 
observational hints [3]. 

On the theoretical side, to generate an observable amount of non-Gaussianity requires 
significant modifications on the standard single-field slow-roll inflation scenario, where all 
the density perturbations come from the quantum fluctuation of the infiaton. Perhaps 
the simplest modification is to introduce an additional light scalar field, since the density 
perturbation of such a scalar field always contains some amount of non-Gaussianity. For 
instance, a curvaton field, which accounts for the curvature perturbations, can generate 
large non-Gaussianity [H [5], [6]. It is also possible to generate large non-Gaussianity with 
negligible contribution to the curvature perturbation, and we call such a scalar as an 
ungaussiton [7]. Moreover, a light scalar may be stable and contributes to some fraction 
of the dark matter of the Universe, as in the case of axion [8]. In this case, the axion 
has a CDM isocurvature perturbation, and as a consequence, the non-Gaussianity of 
iso curvature- type is generated P [TU | [TTl \T2\ . 

Non-Gaussianity from a light scalar, particularly in the curvaton scenario, has been ex- 
tensively studied in the literature [ISHH], but most of them considered a simple quadratic 
potential or a potential slightly modified from a quadratic one [T5t [16] . In this paper we 
study how the resultant non-Gaussianity is affected in a case that the scalar potential is 
largely deviated from the quadratic potential. One example is a pseudo Nambu-Goldstone 
(pNG) boson, and the other is a scalar rolling from the origin in a symmetry breaking 
potential. The former example is particularly interesting, since the flatness of the poten- 
tial is protected by the symmetry, which enables the pNG boson to have large quantum 
fluctuation. Thus a pNG boson is a good candidate for the curvaton [1^. Thus it is 
worth studying the non-Gaussianity arising from the fluctuation of the pNG boson, which 
has a cosine-type potential. Despite its importance, less attention was paid to the non- 
Gaussianity in the pNG curvaton model. In particular, if the pNG boson stays initially 
around the top of the potential during inflation, the predictions on non-Gaussianity are 
expected to change compared with those in the standard curvaton scenario with a simple 
quadratic potential and such a case was not studied in the literature, so far. 

This paper is organized as follows. In Sec. [2] a general formalism for calculating the 
non-Gaussainity in the curvaton scenario is briefly reviewed. In Sec. |3] we study the 
non-Gaussianity from pNG type curvaton scenario. Sec. [His devoted to conclusions. 
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2 Non-linearity Parameter 



According to the 6N formalism [181 [El 120], the curvature perturbation ( on sufficiently 
large scales on the uniform density slicing is evaluated by 

( = N{U,tf-x)-N{U,tf), (1) 

where N{ti,tf]x) is the e-folding number from the initial spatially flat hypersurface at 
t = tj, to the final uniform density slicing at t = tf along the world line x =const., and 
N{ti,tf) is the background e-folding number. We take the initial fiat slicing at a time 
when a scale of interest exits the horizon during inflation, and set the final uniform density 
slicing in the final radiation dominated epoch before the cosmological scales re-enter the 
horizon. ( is known to be a gauge-invariant quantity and time-independent in the absence 
of non-adibatic pressure [21]. Assuming that there are hght scalar fields, {(pa}, during 
inflation, ( can be expanded as 

C = NMa + ^Nab6(l)aS(pb + ■■■, (2) 

where Scpa denotes the fluctuation of 0^ during inflation at t = tj, and A^^ = dN/dcpa, and 
so on. The power spectrum and the bispectrum of ( are written as 

(Ck^Ck,) = {2rcmh + h)Pdh), (3) 
iCk^Ck^Ck,) = {2nmh + h + h)Bc{h,h,h). (4) 
The non-linearity parameter /nl is defined by 

Bcih, k2, h) = ^/nl [Pdh)Pc{k2) + Pc{k2)Pdh) + Pc{k,)Pdki)] . (5) 

In a squeezed configuration such that one of the fcj's (i = 1, 2, 3) is much smaller than the 
other two, /nl can be expressed in terms of Na and Nat as 



where kf, = min{fci, k2, k^}, and we have introduced an infrared cutoff L, which is taken 
to be a few times larger than the present Hubble horizon scale |22l [23]. We have here 
defined 



^s, - ^Psm - ( ^ ) , (7) 



where ifjnf is the Hubble parameter during inflation. 

Let us consider a simple curvaton model where a light scalar field a is responsible for 
the total curvature perturbation, and the inflaton has negligible contributions to the 
bispectrum. Then the expression is simplified to 

5/NL - (iv2 + • 
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Assuming that the oscillating a field behaves like non-relativistic matter, N^, and Nf,a are 
given by [H] 




2R 




2 



3 



(9) 



where a^s denotes the field value at the onset of its oscillation and the prime denotes a 
derivative with respective to the initial position at t = tj, denoted by (Tj, and R roughly 
denotes the fraction of the curvaton energy density at the time of its decay, 



In the case of quadratic potential, o"os = /3crj with an C(l) constant (5 = (2/mtos)^''^r(5/4) 
X Ji/4(mtos) {Jn{x) is the Bessel function) where tos denotes the epoch at which a begins to 
oscillate, and hence a'^^/aos = ^/(^i- If the curvaton potential receives correction from non- 
quadratic terms, the resulting non-Gaussianity may significantly change [15]. Assuming 
that the curvaton potential is quadratic and the curvaton fluctuation is dominated by the 
linear part, i.e., \5a/ai\ ^ 1, we obtain 



where we have neglected the contribution to the curvature perturbation from the inflaton 
fluctuation. Thus a small value of R is necessary to generate large non-Gaussianity. 
Note that the above discussion relied on an assumption that the curvaton behaves like 
non-relativistic matter, and this means that the curvaton potential is (almost) quadratic. 
However, we are interested in a case where the curvaton oscillates in the cosine-type 
potential, and the above formula cannot be applied to such a case. In the next section 
we study the non-Gaussianity in the pNG curvaton scenario. 

3 Non-Gaussianity from pseudo Nambu-Goldstone bo- 
son 

As is well known, a large non-Gaussianity can be obtained in the curvaton scenario, only 
if -R ^ 1. However, as already noted in the introduction, one of the difficulties to realize 
the curvaton scenario is that scalar fields other than the infiaton generically have masses 
of order of the Hubble scale in supergravity and hence their quantum fiuctuations are 
highly suppressed. In order to forbid such Hubble mass terms, one needs some special 
infiation model such as the D-term infiation model [24], or tuning of the coupling between 
the curvaton and inflaton. One of the most promising ways to achieve this is to identify 
the curvaton as a pNG boson associated with spontaneous symmetry breaking, since it 
does not receive a Hubble mass correction. Therefore a pNG boson is a good candidate 
for the curvaton [17] . 




(10) 




(11) 
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Figure 1: Schematic pictures of the quadratic and axion-type potentiaL 



Consider a case that a complex scalar $ has a vacuum expectation value and spon- 
taneously breaks U(l) symmetry. If the symmetry is exact, a NG boson corresponding 
to the angular direction of $, denoted by a, remains massless. But if the symmetry is 
explicitly broken by a small amount, a feels an axion-type potential with a symmetry 
breaking scale /, 



2 r2 



V{a) = m'f 



a 



(12) 



where m is the mass of a field. 

Let us suppose that the symmetry is already broken during inflation, and is not re- 
stored after that. Then a freezes at the initial value during inflation (cij) until it starts to 
oscillate. In principle (Xj can take any value between and 2-7?/. If a happens to sit near 
the origin {\ai\ <^ vr//2), its potential is approximated by 



(13) 



and hence it reproduces the usual results of the curvaton models with a quadratic potential 
(see Fig. [1]). However, a may be relatively close to the top of the potential where the 
above approximation breaks down. This case was not studied in previous works. One 
complexity arises due to a fact that the oscillating a does not behave like non-relativistic 
matter since the potential is significantly deviated from the quadratic one in this case. 
We have performed a numerical calculation to evaluate and N^a- We have solved 
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a set of equations 



dV 

^ + ^Ha = - — , (14) 



da 



H' = ^^iPr + P^), (15) 
p. = + Via). (16) 

from the initial time ti, taken to be well before a begins to oscillate, to the final time tf, 
taken to be the time of a decay, under a sudden decay approximation. For various initial 
conditions ai, we calculate the total e-folding number N{ai), and its first and second 
derivative, N„ and N„„. We calculate /nl using Eq. ([H]), where we assume that ai is 
chosen in such a way that the first term dominates over the second one in Eq. ([8]). As cxj 
approaches the top of the potential for fixed i^inf, the e-folding number N diverges since 
a stays there for longer and longer duration in that limit and the perturbative expansion 
breaks downj^ 

One may expect that the sign of /nl changes around a, = vr//2, reflecting the change 
of the sign of the second derivative of the potential. The naive estimation goes as follows. 
Noting that the curvature perturbation generated by the curvaton a is given by 



^ I5p^ R 
3 ptot 3 



(17) 



where ptot is the total energy density of the Universe at the curvaton decay and the prime 
denotes the derivative with respective to a, we naively expect 

R sin(a,//) ^R cos(a,//) 

- 3 /[I - cos(a,//)] ' 3P[1- cos(a,//)] • ^ ^ 

Thus using Eq. ([H]), we would obtain 

r _ 5 cos{ai/f)[l-cos{ai/f)] 

/nl ^ ■ 2f /r-, • (19) 

2R sm {ai/f) 

In the limit ai <^ /, the result for the quadratic potential /nl — 5/(4_R) is reproduced. 
On the other hand, this expression shows that /nl changes the sign for ai > (7r/2)/. In 
particular, if the curvaton is initially placed near the top of the potential, we would obtain 
/nl — ~(//^<^)^(5/-R), where Act = tt/ — ai. However, note that this expression ignores 
the kinetic energy, and we show that this naive estimation is not correct. 

The results of the numerical calculation are shown in Figs. [2] and [31 In Fig. [2] we have 
plotted N„ and N„„ in Planck units, as a function of the initial position ai. The left panel 
of Fig. [3] shows /nl calculated from Eq. ([S]) assuming that the curvature perturbation is 



^In reality, the e-folding number does not go to the infinity, and instead, domain walls will appear, 
making the estimate on the non-Gaussianity more complicated. 
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Figure 2: N^j (top) and N^fj (bottom) in Planck units as a function of initial position 
of the curvaton, o",, for the quadratic and cosine-type potentials. Here we have set / = 
5 X 10~^Mp. We also show the naive estimation of /nl for the cosine-type potential. 



6 



dominantly generated by the curvaton (A^o- ^ ^<t>)- For comparison, we have also shown 
the numerical result for the quadratic potential f[T^ . Here we have set the symmetry 
breaking scale as / = 5 x 10~^Mp and the decay rate of a as T^r = 10~'^m, and assumed 
that the reheating temperature after inflation is so high that the Universe is dominated 
by the radiation when a starts to oscillate. It is seen that the value of N^r and N^ra 
for the cosine-type potential significantly deviate from the quadratic one, and also from 
the naive estimation for the cosine-type potentially Surprisingly, however, the resulting 
non-Gaussianity for the cosine-type potential is very similar to that of the quadratic one, 
irrespective of the initial field value. This is due to an accidental cancellation between 
and Nf^a in ([8]). Actually, the sign of /nl is always positive even if the initial position of 
the curvaton field is located near the top of the potential. This does not match with the 
previous naive estimation. The resolution lies on an implicit assumption that the curvaton 
begins to oscillate at H = m and afterwards it behaves exactly like non-relativistic matter. 
Both are incorrect in general. If the initial position of the curvaton is near the top of the 
potential, the oscillation epoch is significantly delayed, which causes more contribution to 
the e-folding number, leading to the enhancement in iVo- and A^o-o- (see Fig. [2]). Therefore 
the positive /nl is obtained even in this case. This should be contrasted to the case 
studied in Ref. [U]. They studied the curvaton potential with a small positive correction 
to the quadratic potential, and obtained a negative value of /nl- This can be understood 
as a result of the suppression of the growth of the e-folding number due to the positive 
correction termjf] 

In the above arguments, we have assumed that the curvature perturbation is domi- 
nantly generated by the pNG boson. The situation changes if the inflaton accounts for the 
curvature perturbation. In the right panel of Fig. [3] we show the non-linearity parameter 
/nl in the presence of the contribution of the inflaton to the curvature perturbation, set- 
ting A"^ = 1 in Planck units. In this case, the non-Gaussianity in the cosine-type potential 
can be orders of magnitude larger than that in the quadratic potential. 

This is a good news for those who want to produce large non-Gaussianity in the 
curvaton scenario, since one needs not to tune the initial value of the curvaton field. Once 
there are pNG fields during inflation, it always predicts positive value of /nl irrespective 
of its initial position. Of course, in order to obtain observable level of non-Gaussianity, R 
should be smaller than about 0.1, and it requires a tuning of the decay rate of a. However, 
in the pNG scenario both the amplitude and the decay rate of the curvaton can be in 
principle determined by an underlying physics, not by rather arbitrary initial conditions. 
Actually the decay rate of the curvaton is expected to be ~ rn? / p and its initial 

^ In order to reproduce the observed density perturbation, the inflationary scale should take an 
appropriate value for each point. For example, N^j < 1 requires unacceptably large Hubble scale during 
inflation, Hmi ^ 10'^'' GeV. However, this is just an artifact of our choice of for our numerical 
calculation. The inflationary scale can be sufficiently small if we take a smaller value of To-. 

^ It was noticed in Ref. [16] that a negative correction to the quadratic potential leads to an enhance- 
ment on /nl, in the case that the correction term is small enough and can be treated as a perturbation. 
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Figure 3: Tlie value of non-linerity parameter /nl as a function of initial position of the 
curvaton, cxj, for the quadratic and cosine-type potentials, for N„ ^ (left) and A^^ = 1 
(right). Here we have set / = 5 x lO^^Mp and To- = 10~^m. We also show the naive 
estimation of /nl for the cosine-type potential. The value of /nl decreases (increases) 
as Fg- becomes smaller for the left (right) panel, although the relation between the three 
lines will not change. 



amplitude is ~ /. Thus R can be estimated as 



«-0.3|^V / YYi^r', (20) 



105 GeVy VlO^^ GeVy \ m J 

where is the reheating temperature after inflation. Although it may be difficult to 
distinguish such a scenario from standard curvaton models, it is important to keep in 
mind that the detection of non-Gaussianity may be closely related to an underlying high- 
energy physics. 



4 Discussion and Conclusions 

In this letter we have investigated the non-Gaussianity from a light scalar with a potential 
largely deviated from the quadratic one, focusing on the pNG curvaton scenario, motivated 
by the fact that the required flatness of the potential is naturally ensured. What we have 
found is that if the potential significantly deviates from the quadratic one, naive estimation 
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on the non-linearity parameter /nl is not correct in general. This is because the timing 
when the oscillation begins and the equation of state during the coherent oscillations are 
different from those in the quadratic potential. Interestingly, if the curvaton accounts for 
the observed curvature perturbations, the non-linearity parameter is quite similar to that 
of the usual quadratic one, irrespective of the initial position of the curvaton field. If the 
inflaton is responsible for the curvature perturbation, the non-Gaussianity is enhanced 
in the pNG curvaton case. Therefore, the pNG curvaton scenario is a good example for 
realizing a large non-Gaussianity. 

Although we have focused on the pNG potential, a similar analysis can be applied to 
other hilltop potentials. For example, consider the potential given by 

Via) = - (21) 

with self-coupling constant A(> 0), which has an extremum aX a = m /TA. If the a field 
rolls down toward the origin from near this extremum, the resultant non-Gaussianity is 
similar to that predicted in the corresponding quadratic potential, as in the pNG case. 
Therefore our results indicate that if the curvaton has a similar hilltop potential which 
changes to the quadratic one at the potential minimum, the non-Gaussianity can be 
estimated approximately by the corresponding quadratic potential, if a is dominantly 
responsible for the total curvature perturbation. 

Finally we comment on the non-Gaussianity generated by the axion. In Ref. [Hj, non- 
Gaussianity from the axion was studied focusing on the case that the axion initially stays 
near the potential minimum and the potential is well approximated by the quadratic one. 
According to the result of this paper, if the axion is placed near the top of the potential, 
the resultant non-Gaussianity as well as its isocurvature perturbation are enhanced. 
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